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Abstract 

The thermodynamic Bethe Ansatz equations arising in the context of the AdS5/CFT4 
correspondence exhibit an important difference with respect to their analogues in relativistic 
integrable quantum field theories: their solutions (the Y functions) are not meromorphic func- 
tions of the rapidity, but live on a complicated Riemann surface with an infinite number of 
branch points and therefore enjoy a new kind of extended Y-system. In this paper we review 
the analytic properties of the TBA solutions, and present new information coming from their 
numerical study. An identity allowing to simplify the equations and the numerical implemen- 
tation is presented, together with various plots which highlight the analytic structure of the Y 
functions. 



§ 1. Introduction 

In 1997 Maldacena [1] (see also [2]) conjectured a correspondence between grav- 
ity and conformal gauge tlieory of strong/ wealc type and, in particular, the duality 
AdS5/CFT4 between the supersymmetric string theory of type IIB on a curved space 
AdSs X and TVT = 4 super Yang-Mills (SYM). 

A turning point for the study of the correspondence AdSs / CFT4 was the discovery 
of integrability in the 't Hooft limit, both for the string theory [3l S] and the conformal 
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field theory CFT4 |5]. The 't Hooft planar limit is defined by the scaling of the colour 
number N ^ 00 and the SYM coupling qym while keeping the coupling Nqym = 
A = 47r^(7^ finite, so that g is proportional to the string tensioili]. In this limit only 
planar Feynman diagrams survive [6]. 

The Bethe Ansatz (BA) method was originally introduced by H. Bethe [7j in 1931 
to study a one-dimensional quantum system of spin 1/2 and determine its eigenvalues 
and eigenstates. This system is also known as the Heisenberg XXX model and, since 
the scattering of the possible excitations can always be factored in terms of two-body 
amplitudes, it has the feature of being integrable. 

In the = 4 SYM context, a set of Bethe Ansatz-like equations were formulated 
by Beisert and Staudacher [8] (see also [9]) as a useful tool to compute the anomalous 
dimensions of SYM composite single trace operators or the energy of the quantum 
states of the string. However, the equations emerging in this framework do not take 
into account the wrapping effects due to the finite number of elementary operators in the 
trace [TOl [11]. Therefore, the Beisert-Staudacher's equations lead to exact results only 
for asymptotically long operators and are the analogue of the Bethe- Yang equations 
for scattering models describing the quantisation of momenta for a finite number of 
interacting particles living on a ring with very large circumference \12\ [T3] . 

A first important step towards the partial solution of the wrapping problem was 
made in [1^ . Adapting the formulas for the finite-size Liischer corrections plSj [T6] to the 
AdS5/CFT4 context, Bajnok and Janik were able to predict the four loop contribution to 
the Konishi operator. The result was readily confirmed by the perturbative calculations 
of [17]. The ideas proposed in [H] can be applied to higher perturbative orders and 
twists p!8l [To] but technical complications increase dramatically with the loop order 
and, through this approach, exact all-loop results are probably out of reach. 

An alternative tool to find the finite-size corrections to the ground state energy 
in 2D integrable models was proposed many years ago by Aliosha Zamolodchikov [20] 
modifying the thermodynamic Bethe Ansatz (TBA) whose original formulation traces 
back to Yang and Yang [21]. The TBA method leads to a set of coupled non-linear 
integral equations governing, as the parameters of the model are changed, the ground 
state energy of the theory on a cylinder exactly and non-perturbatively. 

As discussed in [22], consider a two-dimensional euclidean quantum field theory 
defined on a torus with circumferences R and L, as shown in Figure [TJ The quantiza- 
tion of this theory is equivalently obtained by considering two alternative hamiltonian 
schemes. In the direct scheme the system is quantised along the direction L and the 
hamiltonian "Hl propagates the states in the time direction R. In the mirror scheme 
the system is described by the hamiltonian T-Ln and the roles of L and R are swapped. 



Another definition for tfie coupling g is also widely used, so that in many works is found the relation 
A= IBttV- 
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Figure 1. The two hamiltonian descriptions of a 2D euclidean model. 

The two descriptions lead to the same partition function 

(1.1) Z{R,L) = Tr(e--^^^) = Tr(e--^^«). 

Sending — )■ 00, in the direct description of the theory the partition function is domi- 
nated by the smaUest eigenvalue Eq{L) of Hl- 

(1.2) Z{R,L) = Tr{e-^^^) ~ e'-^^o^^^ 

While, in the mirror description, the R 00 limit corresponds to the thermodynamics 
of a one-dimensional quantum system with hamiltonian T-Lfi defined on a volume R at 
temperature T = In this picture, the partition function becomes 

(1.3) Z{R, L) = Tr(e-^^«) ~ e'^^f^^l 

where f{L) is the free energy per unit length at equilibrium. Comparing ( II. 2p with ( II. 3p 
we see that the ground state energy of the direct theory is related to the equilibrium 
free energy of the mirror theory: 

(1.4) Eo{L) = Lf{L). 

Starting from the Bethe-Yang equations for the mirror theory, f{L) can be found [20] 
adapting the method proposed in [21] for a one-dimensional system of bosons with 
repulsive delta-function interaction. As a result of the TBA procedure, the exact finite- 
size ground state energy Eq[L) is written in terms of the pseudoenergies £„, solutions to a 
set of coupled non-linear integral equations. In [231 1211 IISI |26] the method was extended 
to excited states and an alternative but equivalent approach was proposed in [271 [2HI [29| 
(This alternative tool has also been used in the AdS/CFT setup in |30l [311 [321 [33]). 
Starting from the mirror version of the Beisert-Staudacher's equations [34| 155] , the 
ground state TBA equations were independently derived in [351 ISZl [SH] and precise 
conjectures for particular excited state sectors were made in [321 [23 [ID] • 
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Because the equations for excited states are expected to be the analytic continuation 
of the ground state equations [211 125] , the understanding of the analytic structure of the 
TBA system is essential in the attempt to establish the results of [39l [371 HD] rigorously 
and to generalize them to other sectors of the theory 

The rest of this paper is organised as follows. Section [2] contains the TBA equations 
of [361 EH] . The pseudoenergies have an infinite number of square root branch cuts 
and, basically, the AdS5/CFT4 TBA equations are equivalent to a set of functional 
relations containing also information on the discontinuities across the cuts [41J: the 
extended Y- systems described in Section [31 Following [H], Section [H discusses briefly the 
interpretation of the TBA as dispersion relations and, using byproduct identities of [H] , 
a variant of the TBA equations more suitable for numerical integration is proposed 
in Section [31 Sections [6}[8l contain our preliminary numerical results for the ground 
state energy and a study of the pseudoenergies in the complex rapidity plane. Finally, 
Section [9l contains our conclusions and the S-matrix elements corresponding to the TBA 
kernels are listed in Appendix \K[ 

§ 2. The TBA equations 

The TBA equations with arbitrary chemical potentials {Ha} are [351 1221 [3H] : 

SQiu) = fiQ+LEqiu) - ^Lq/ * (Pq'q{u) 

Q' 



(2.1) 



M 



(2.2) 



Q M 



(2.3) = ^(HK) +5Z^HM) +4"^*7c ^k{u), 

M 



(2.4) = - * + E(^(:|M) - ^Sm)) * Mu). 

Q M 

where L G N is the inverse of the temperature, while the sums are taken as ^q.^^ 2' 
Ek=1' with 

Eq{u) = In ^ ^, x{u) = ^ - ^1 - T ' < 0' 

x{u + If) \ 2 V 4 / 
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^=^^ r 

^ ( y — ^ ^ y \- 



2 < 



Figure 2. The contour 70. Figure 3. The contour 7x 



(2.5) Ya{u) = e^'^(^), La{u) = ln(l + 

and the symbols and '*-y' denote the convolutions 



J" * Mu) 



I dz J^{z) (t){z,u), T (p{u) = (b dz ^{z) (p{z,u) 

Jr J-y 



The kernels are derived from the matrix elements Sabi^iu) listed in Appendix Rl using 
the relation: (j)ab{z,u) = i^^-j^lnSabiz^u). (Pq/q in ( 12. ip can be written as: 

(2.6) (Pq'q{z, u) = -(pQ'Qiz - w) - 2K^,q{z, u), 

with [36l[371|38lll2lll3lllllll5lllT] 

2.7 (Pq'qM = -InSn'oM, Kp(z-t) = ^ In (v4, 

(2.8) K^,q{z,u) = ^.^\nT.Q,Q{z,u) = ds (I)q, ^y{z, s) dt K^is -t)(/}y^Q{t,u). 



2711 du 
1 d 



The contours of integration % and 7x are represented in Figure [2] and Figure [3l respec- 
tively. The functions Ya{u) = e^°-^'^^ solutions of equations f l2.1tl2.4p live on branched 
coverings of the complex u-plane with an infinite number of square root singularities 
u = {±2 + im/g} with m G Z. For the mirror AdS5/CFT4 theory under consideration 
-on the sheet containing the physical values of the Ys- all the cuts are conventionally 
set parallel to the real axis and external to the strip |Re(w)| < 2. This will be referred 
to as the reference or first sheet. Table [U shows the location of the square branch points 
for the various Ys. Finally, we shall denote with Y^'^^_^{u) (or simply Yy'^\u)) the first 



Function 


Singularity positions 




u = ±2 + i^, J = 0,±1,±2, ... 


y^"^ (u) 


u = ±2 + i^, J = ±M,±{M + 2),±{M + 4),... 




YMiu) 



Table 1. Square root branch points for the Y functions. 
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sheet determination of Yy and with Yl^y^'_^^{u) = Yy (u^,) its second sheet evaluation 
obtained by analyticaUy continuing u to through the branch cut u G (— c>o, —2) (see, 
Figure E]) . 



§ 3. The extended Y-system 

The Y-system |l6l IHl SH SH] for AdS5/CFT4 was conjectured in ^J, rigorously 
derived in [36l[37l|38] and it is associated to the diagram represented in Figure HI Setting 

(3-1) 

M M 



with Ca// AT = 25m,n — Amn-: ^i,m = 82,m, Anm = 5m,n+i+Sm,n-i and Amn = Anm, 
the Y-system with arbitrary chemical potentials is: 



5q.i-1 



(3.2) YQ{u-l)YQ{u + l)=AQl[il + YQ,{u)f^<^'l[ ^ ^ 

Q' « M _L 1 



(3.3) y; 

(3.4) 



A 



(.) v^ + ^(S)H 



1 + 



{w\My 



(3.5) 



AT 



1 + 



1 + 



Y, 



(al + ) 



U — 



^*-(«|M) 



Hat ( 1 + -^(llAT) 



1 + 



M + l 



(u) 



l + Y, 



(a) 



iy\+) 



[u) 



In relativistic integrable models the Y functions are in general meromorphic in the 
rapidity u with zeroes and poles both linked to 1 + Ya zeroes through the Y-system. 
Unfortunately, the situation for AdS5/CFT4 is further complicated by the presence of 
square root branch cuts inside and at the border of the fundamental strip |Im(u)| < 
1/g. According to the known Y — )■ TBA transformation procedures, together with the 
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(w|N) 



a=l 



a=2 




(w|N) 



Figure 4. The Y-system diagram corresponding to the AdS5/CFT4 TBA equations. 



asymptotics of the Y functions this extra information on their analytic behaviour should 
be independently supplied. 

However, the crucial discontinuity information is stored into functions which depend 
non-locally on the TBA pseudoenergies [3Sl IH] and thus on the particular excited state 
under consideration. 

The main objective of the paper was to show how this problem can be over- 
come and that all the necessary analytic information can compactly be encoded in the 
basic Y-system ( 13. 2113. 5 p extended by the following set of local and state-independent 
discontinuity relations. Setting 



(3.6) 



A{u) = [InYiiu] 



1+1 ' 



then A is the function introduced in [38] and the discontinuity relations are: 
(3.7) 



[A] 



±2N 



0=1,2 



In 1 



1 \ 



Yr 



(y|T) 



N 

±2Ar M=l 



In 1+ 



1 



Yr 



(a) 



(v\M) 



±(2N-M) 



(3.8) 



'Y, 



(a) 



In 



iy\-) \ 



±2N 



N 

-E 

Q=l 



In 1-f 



Y, 



Q. 



±{2N-Q) 
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with iV = 1,2, 



oo and 



(3.9) 



In 



±1 



1 _L 

^ + ^'^(y\-) 



InK 



(a) 



±1 



In 



(a) 



(J/I-) 



Vi 4- 



where the symbol [f]z with Z G Z denotes the discontinuity of f{z) 
(3.10) [/]^ = hm f{u + iZ/g + ie) - f{u + iZ/g - ie), 

on the semi-infinite segments described hj z = u + iZ/g with u G (— oo, —2) U (2, oo) 
and the function [f{u)]z is the analytic extension of the discontinuity ( I3.10p to generic 
complex values of u. Finally, A(u) has an additional constant discontinuity running 
along the imaginary axis: 



(3.11) 



A{iv + e) - A{iv - e) = z2L7r, {v G 



In conclusion, while the Y functions are defined on an infinite sheeted Riemann surface, 
the basic Y-system ( 13. 2113. 5 p connects points lying on a single reference sheet missing a 
huge amount of analyticity data. To recover this information it is necessary to extend 
( I3.2tl3.5p by including relations among different Riemann sheets. This is precisely the 
information contained in discontinuity relations and, to highlight this fact, one could 
write equations ( I3.6ti3.9p in a more explicit functional form involving points from differ- 
ent sheets. For example (cf ( I3.9P ): 



(3.12) 



y;7,,(w, ±^) Vi + i/y/"l, 

('u;|l)V * g) ' / (y|+) 



{v\l)\ g> 



y/"\(w, ± ^) 

{v\l)\ * g> 



{y\-) 



(u) 



where is the second sheet image of u reached by analytic continuation through the 
branch cut u G (— oo, —2) (see. Figure E]). 




Figure 5. The second sheet image of u. 



4. TBA equations as dispersion relations 



As already mentioned in the introductory section, contrary to the more studied 
relativistic invariant cases, the transformation from Y-system to TBA equations is by 
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no means straightforward since the local form of the AdSs / CFT4-related Y-system (13. 21 - 
I3.5p alone does not contain information on the branch points and it is almost totally 
insensitive to the precise form of the dressing factor Eg/g defined through equation (12 .Sp . 
The purpose of this section is to give some hints on why the extension of the basic Y- 
system by the discontinuity relations (I3.7tl3.9p resolves this problem completely. The 
interested reader is addressed to [H] for a more complete discussion on this important 
issue. Consider the following equation directly descending from (12. 4p : 

(4.1) l^iylyt)i^)/ylyt)(^)) = - Y. I dzLQiz){Kiu-iQ/g,z)~Kiu + iQ/g,z)), 



Q 



where 



(4.2) K{z,u) ^ 



27riV4 - z'^ z - u 

Then, it can be shown that the functional relation (13. 8p on the discontinuities 

N 

(y|-)/ 



(4.3) h(^(';-)/4il))l.,„ = -EM± 



±2N — ' ^ ^' 

Q=l 



combined with some more general analyticity information, is equivalent to (14. ip . This 
result comes from the observation that the quantity: 



is analytic at the points u = ±2, but it still has an infinite set of branch points at 
u = ±2 ± i2N / g with G N. From the Cauchy's integral theorem we can first write 

^^^^ in(i',';i,(")/K<;t,w) _ ^ m (r,';i,(.)/F,'„|i,(.) 



V4 - 27rz {z - w)V4 - z"^ 

where 7 is a positive oriented contour running inside the strip |Im(M)| < 1/g, and 
then deform 7 into the homotopically equivalent contour Tq represented in Figure [6] 
as the union of an infinite number of rectangular contours lying between the branch 
cuts of (14. ip . Restricting the analysis to the solutions that behave asymptotically as 
ln(Yj.^J|^-j/Yj.^J|^P — )■ 0(1) uniformly as \u\ — )■ 00 the sum of the vertical segment con- 
tributions vanishes as the horizontal size of the rectangles tends to infinity. Then, 
equation (14. ip is recovered after inserting relation (14. 3 p and observing that several can- 
celations take place between adjacent integral contributions. In [H] the entire set of 
TBA equations was retrieved starting from the extended Y-system and adopting cer- 
tain minimality assumptions for the number of logarithmic singularities in the reference 
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Figure 6. The deformed contour To- 



Riemann sheet. These are precisely the analytic conditions fulfilled by the ground state 
TBA solutions as fXy ^ in. Here, we would like to mention that is possible to show [50] 
that these assumptions can be relaxed, and that logarithmic singularities of the and 
La functions can lie on the first sheet without affecting the derivation provided they are 
far enough from the real axis and organized in complexes as required by the Y-system 
relations. We can summarize the situation at zero chemical potentials in the following 
proposition. 



Proposition 4.1. 

Let the singularities on the first sheet be organized in complexes such as the follow- 



ing: 



(4.5) 



or 



Yq{u\ 



{v\Q-i)y^i + 



9 



-1; 

0, 



Y, 



I M+l) ((""■■ 

(4.6) 



(-1) 



(a) 



(v\M-l) 



(-1) 



Y, 



(a) 



2i. 
9 



^{v\M)y"'2 > - 

0; Ym+i{u2 H — ) = oo^ ( double pole ), 
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or 



(w\M)('^3 



Y, 



(a) 



{w\M+l) 



(-1) 



+ 



(w\M)\ 3 ' g . 



Y, 



{w\M-l) 



(-1) 



+ 



-1; 
= 0, 



(4.7) 



or 



+ -) = 0; Yiiuf'^ + -) = oo'' ( double pole ) 



,(-1) 



^(y|-)l«4 . 

(-1) 



Yr 



(y|-)^"4 



2k 



where, by definition, in equation \4-5^ - ^{v\o) ~ ^{y\-)' ~ ^' ^^'^ equations ( [^.6' 
IS- ^HO) = l/^(Ho) = ^(y|±)- 



Moreover, let the zeroes u'^^^ = g outside the strip \Im{u)\ < 

r/ien t/ie extended Y-system implies the ground state TBA equations. 



The ground state equations are modified only when the two points u^~^^ and u^~^^ + ^ 
lie on different sides with respect to the real axis, which can be interpreted as the result 
of a singularity having crossed the integration contour. 

As we shall see in the following sections, the ground state TBA solution at /x^ = 
shows indeed several of these complexes of singularities, in the reference sheet but far 
from the real axis. 



§5. A useful variant of the TBA equations 



To facilitate the numerical implementation of the equations, the following auxiliary 
functions E^"\ G^'^\ T^"^ and U^""^ are introduced: 

= e<j'_,(«)-.!^,V,(u) 

(5.1) 

Q 



^ dz Lq{z) [k{z -i^,u)-K{z + i^,u] 
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(5.2) 



M 

- V / dzLQ{z)(t)Q{z-u), 



(a) 

y 



(5.3) T^-\u) = L\;l^iu) - L\;l^iu\ u^-Hu) = l[;1^{u) + l[;1){u). 

Using these definitions the TBA equations ( I2.m2.4l) can be recast in the form: 
(5.4) 

4l\K)i^^=^'MK)+Yl [ dzL\i^^,^iz)cPMKiz-u)+ [ dzT^-\z)4>K{z-u), 

(5.5) 

= ^(^IK) + E / ^^^("|M)(^) <I^mk{z -u)+ [ dz {T^-\z) + E^-){z)) M 
-Y] / dzLQ{z)tpQK{z-u), 



z ~ u] 



Q 



(5.6) 



^q(w) = +^^q(^^) / ^^^Q'(^)^Q'Q(^~^) 

+ J2rdz i^^^\z) + ^f/(")(z)) (k(^, ^ - ^f ) - i^(;2, u + i^ 

2 

-^E / dzT^"Hz)Mz-u) + J2 [ dzL[':l^^{z)ijMQiz-u) 



with 



(5.7) XQiz,u)=Q dt + J^ dt^ Kf\z-t)(K{t,u-i^)-K{t,u + i^)y 

where the last term of (15. 6|) is obtained from the last term of ( 12. 6p using the foUowing 
important relation discussed in Appendix D of 



Y,Lq'*K^,q{v) = Y,Lq'* I ds(t^^^)^y{s) I dtKP{s-t^y,Q{t,v) 

O' Q' "'Tx -'Tx 

-I ds lnyj")(s) / KP{s-t)(Py^Q{t,v). 
J 1^ J 



(5.8) ^' ^ 
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To further optimize the implementation of the equations, we have minimized the number 
of terms that depend on the two variables z and u separately H. The last term appearing 
on the rhs of (15 .Sp was rewritten considering relation (E8) in |41] : 



(5.9) 



-V / dzLQ{z)(\)Q^(^^\K){z,u) = / dzE 

-V / dz LQ{z)tpQKiz - u) 



with 



M-1 

(5.10) tpQM{u) = (pQ-M+2j\ 

J=0 



( step 2\Q + M\-2 

E j=\Q-M\ (t>j{u) ] Q> M; 



u) = < 



stop 2\Q + M\-2 

. E j=\Q-M\+2 <^jN ; Q<m. 



Similarly for the penultimate term on the rhs of (15. 6p we have used 
(5.11) 

„, i\/r JR ^ J —2 



Z- 

9 



a,M 



Z - u), 



with 



(5.12) ^i-)(u) = Y, I dzL^^^^^{z)Mz-u). 

N ^R 

The equations for the pseudoenergies used in the numerical work to be described shortly 
correspond to an appropriately discretized version of (15. 4115. 121) where the infinite sums 
are truncated to a finite number of terms N^ax and the integrals performed only in the 
range {—Zmax',Zmax)- The chemical potentials are set to zero but for those associated 
to the fermionic particles (?/|±) where fiy — /j,^^ = —fJ.^^- The results obtained and 
partially discussed in this paper correspond to the parameters in the ranges reported in 
Table U 



This was particularly useful to optimize the occupation of the computer internal memory. 
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Convergence threshold for La 


e = 1 X 10-^4 


Couphng 


(0,1] 


Scale 


L e (0, 10] 


Integration range {- Zmax] Zmax) 


Zmax 50 


Number of intervals in (—2; 2) 


Np e (1,150] 


Truncation 


Nmax e [2, 15] 


Chemical potentials for (y ±) 


€{(-10,10) x[0,7r)} 



Table 2. Parameter ranges. 



§ 6. TBA numerical solution: the real axis 

The functions La{u), solutions of the thermodynamic Bethe Ansatz equations, are 
represented in Figure [7] to Figure [T21 

Re L(u) 

0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 


-15 -10 -5 5 10 15 

u 

Figure 7. Solutions Re(LQ): L = A, g = 0.5, Hy = 0. 

Figure [7] shows the real parts of Lq{u) for Q = 1 and Q = 2dXL = A, g = 0.5 and 
= 0. The solutions are well localized around the origin, tend exponentially to the 
asymptotic value Lg = and in addition Li{u) ^ L2{u). These are general features of 
the solutions of this TBA: at moderate values of the coupling constant g the difference 
between La{u) and its asymptotic value is always well localized about the origin 
and a given component Lain) strongly dominates over the next, ie La{u) ^ La+i{u). 
This justifies both the relatively small range of integration and the truncation of the 
sums to a finite number of terms Nmax- Unfortunately, both the localization and the 
subdominancy features get -for generic values of Hy and L- worse and worse as g is 











Q=1 
Q=2 
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increased above 1 and in the strong coupling regime one has to push the computer 
internal resources to the extreme. Figure [8] shows the results for the excitations of type 
V and w with M = 1, they also reach the asymptotic values = L^w\i) ~ I 

exponentially fast in the large \u\ region. Figure [9] shows instead the real parts of L{u) 
for the excitations v,w with M — 2, they are smaller and reach the asymptotic value 
L^y\2) — -^(u;|2) ~ i fester that the M — 1 cases. Again we see that the M components 
dominate over the M + 1 ones. 

Re L(u) 
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0.2878 
0.28775 

0.2877 
0.28765 

0.2876 
0.28755 

0.2875 
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u 

Figure 8. Solutions Re(L(„[i)), Re(L(u,|i)): L = i, g = 0.5, fiy = 0. 
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Figure 9. Solutions Re(L(„[2)), Re{L(^^\2))- L = A, g = 0.5, = 0. 
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The imaginary parts of L{u) for the excitations Q, v and w vanish within our 
computer working precision. 



Re L(u) 
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Figure 10. Solutions Re(L(j^|-|-)): L = 4, g = 0.5, = 0. 



Figure [TOl compares the real parts for the (y|±) fermionic particles. These solutions 
have a physical interpretation in terms of particle and hole densities only in the rapidity 
range u G (-2, 2) [36l [371 [38] . 

Indeed, one can see cusps at w = ±2 highlighted in Figure [TTJ Outside the physical 
interval (—2,2) both functions tend to the corresponding asymptotic value In 2. 

The imaginary parts are depicted in Figure [T2] and they vanish for u G (—2,2), 
confirming that only in such interval the pseudoenergies (^) are physically mean- 

ingful. 
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Re L(u) 
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Figure 11. The square root branch point of Re(Ly): L = A, g = 0.5, pLy = 0. 




Figure 12. Solutions Im(L(j^|-|-)): L = 4, g = 0.5, fiy = 0. 



The solutions of the TBA equations expressed in terms of the functions La{u) have, 
in principle, a maximum precision of the order of the threshold of convergence shown in 
Table [2l In fact, the solutions are mainly affected by errors due to the discretisation of 
the domain of integration and the truncation of the infinite sums. The TBA equations 
are non-linear, this makes it difficult to assess a priori the magnitude of those errors but 
we have estimated the precision by comparing solutions at various values of Nmax, Zmax 
and Az. As a result of this analysis, in the range of parameters reported in Table [2], 
the estimated precision is between 10~^ and 10~^. 
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7. E = protected state and scale quantisation 



The ground state energy of the model is obtained from the pseudoenergies for the 
Q-particles as: 



(7.1) 

where 
(7.2) 



Eo{L) 



^ I dudp^ 



Q=l 



27r du 



pQ(u) = gx{u - iQ/g) - gx{u + iQ/g) + zQ, 



is the mirror theory momentum. In A/" = 4 SYM the true vacuum state is protected by 
supersymmetry and from the TBA this state is found by sending the chemical potential 
for the fermionic particles to 



(7.3) 



Eo{L,Hy,g) ^0 \/g,L, 



while all the other chemical potentials are kept at zero. 

Figure [13] shows the occurrence of this condition at different values of the parame- 
ters while Figure [T3] compares the numerical result with the asymptotic estimate given 
in EH: 



(7.4) ~ - 5^ 

Q=l 

with h = —ifiy — TT. 



^^lnfl + 16Q2sin^V^^'^(-) 
2n du \ ^ 2 



general h, large L. 



05 

DC 




Figure 13. Eq versus fiy and g: L = 4. 
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Figure 14. Eq, numerical versus the exact asymptotic result: L = 4, g = 0.5. 



Further, from the point of view of A/" = 4 super Yang-Mills the parameter L should 
be restricted to strictly positive integer values. In fact, L is related to the number 
of elementary operators building the composite trace operator of which we would like 
to determine the anomalous dimension. On the contrary, within the TBA setup L is 
simply the inverse of the temperature and does not necessarily need be quantised. In 
[51] and more generally in [IT] a drastic simplification of the analytic properties of the 
solutions was observed at integer L but, up to now, these mathematical facts have not 
been neatly linked to the gauge theory origin of the equations. 

Indeed, the pseudoenergies change smoothly on the real axis with the scale and the 
chemical potential fiy giving a smooth variation of Eo{L, fiy), as shown in Figures [15] 
and [TBI 

Figure [T^ corresponds to a truncation at Nmax = 2, and we see that the energy tends 
uniformly to zero as fXy ^ in. Notice also that, regardless of fiy, there is a particular 
value of L ^ 2.3 where the curves simultaneously intersect the horizontal axis. 
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Figure 15. Eo versus L, ^y-. g = 0.5, Nmax = 2. 



Figure Uni corresponds to a truncation Nmax = 5, again the ground state energy 
tends uniformly to zero as — t- in and at L ^ 2.07 the various curves intersect on the 
horizontal axis. Increasing Nmax we observed that the intersection point approaches 
L = 2, but unfortunately the program fails to converge at L = 2 for Nmax > 8. 



DC 









Im (i„=.62 

lmn/l.23 








lmn^-1.85 

Imn^=2.46 








Im (iy=3.08 









































































1 2 3 4 5 6 

L 



Figure 16. Eo versus L, ^y-. g = 0.5, Nmax = 5. 
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§ 8. Mapping the complex plane 

In this section we shall describe qualitatively the analytic properties of the TBA 
solutions in the complex rapidity plane. Following [2^ [25] the solutions for complex 
values of the rapidity were first obtained in the fundamental strip |Im(u)| < ^ starting 
from the La{u) computed on the real axis and using the TBA equations as integral 
representations for the pseudoenergies. 

Then, the solutions were calculated in the other strips, defined as 

Strip(+^) = {u I Re(w) G (-00, +00) , Im(w) G (+^, +^)}, 
f8 1") y y 

Strip(-^) = {u I Re(r.) G (-00, +00) , Im(r.) G (-^, -f )}, 

with k = 1,2, ■ ■ - , using the basic Y-system ( I3.2tl3.5p and connecting two points in a 
given strip /c to a point outside it, as exemplified in Figure [T71 In this way the whole 
complex plane can be mapped. Notice that the basic Y-system does not contain a 



strip (+] 
Strip (0) 


L) 


< 


i 

► u+2i/g 


L //T7(U) 

< 
< 


■ +2i/g 

» u 

- +i/g 

( u-i/g 

Re(u) 


Strip (-] 


-) 


< 

r 


> u+i/g 

► u 
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► 

> u-2i/g 

i/g 

2i/g 



Figure 17. Mapping of the complex plane of rapidity u. 
functional equation for , to overcome this problem we used equation (]5.2p 

(8.2) e'^^^^p)=G^-\u)-e^^^_^{u), 
with 

(8.3) ^^"^(^) = E^A? * + 2^"^ 

M 

and £^''(w) = 2L|"|^^^j(u,) — 2L|^|^^^(w,) — Lm{u) and analytically continued it to a 
generic strip in the complex plane: 

(8.4) u G Strip^i'^) ; G^"\u) = J2^M* <Pm{u) + 2^") + 4"^^ T ^n) 

M n=l 
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The last term on the rhs of ( 18 ■4p comes from the residues of the simple pole singularities 
in (Pm{u)- In order to unveil the analytical properties of the Y functions in the complex 
rapidity plane we have introduced the functions 



(8.5) 



1 + + Y^\u)\ 



E [0, 1] 



and considered the special or critical points u^^\u^ ^\u^°°^ such that 



^.6) 



YAiu^'^^) = ; Ya{u 



(-l)^ 



-1 ; YAiu^"^^) oo. 



Table [3] summarizes the values assumed by Fa^\u) in correspondence to these points. 
The functions -F+ emphasize better the and points, on the contrary the 

function F~ are best suited to highlight the points u^^^ and u^~^\ 



u 








Ya{u) 





-1 


oo 




1 

2 





1 


F-iu) 


1 





1 

2 



Table 3. values corresponding to u 



Figure fTSl shows a three dimensional plot of F^ (u) where the square root branch 
cuts emerge and we also see peaks at zero {u^~^^) and peaks at one (u^^^). A picture of 
this kind is certainly nice looking but it is not very easy to interpret. 
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Figure 18. 3D plot of F^{u): L = A, g = 0.5, Hy = 0. 

We therefore consider two dimensional plots with contour lines, in which the value 
of the functions F^"^^ is encoded using a scale of colours. Blue corresponds to F^"^^ = 
while red to F^^^ = 1. Figure fT9l plots -Fj~, where the critical points of type w*-'^^ and 
M*-"^^ and the branch cuts are easily recognisable. The 6aszc Y-system relates poles and 
zeroes between adjacent strips on the reference sheet, so that complexes or strings of 
critical points emerge. As an example of these phenomena, we shall highlight the link 
between points u^^\ u^~^^ and u^°°^ induced by equation (13. 2 p at Q = 1: 

/ y^'^ (u) 

(8.7) Y^{u - ^)Y,{u + ^) = (1 + Y^{u)) ^ ^ 

where we have used the symmetry = ^{y\-) A*l/"^ = 0- Equation (18. 7p shows 

that there are relations between the zeroes and the poles of the functions 

(8.8) Fi, ^2, (1 + Y2), y/„V-v (1 + ^/'^ ' 



(y|-)' ^ '{y\-)' 
(2/1-)' 



From Figure[2U]we see that the functions Ft^,-. {u) vanish at w = w*- ~ (3.7, 4.4) G 
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Re(u) 



Figure 19. Contour plot of (u): L = 4, g = 0.5, fiy = 0. 

Strip (in blue). At this point the rhs of equation ( 18. 7p tends to infinity, this im- 
phes that F^{u) = 1 at m = w^"") ^ (3.7,6.4) e Strip(+^) (a double pole of Yi) which 
can be clearly spotted in Figure [21] (in red). Consider now Figure [211 the functions 
F^{u) show the existence of zeroes at u = u^~^'> ^ (3.1,6.5) G Strip*-"*"^^ (in blue) so 
at that point the Ihs of (j8.7l) is equal to one. In turn, this implies that F2 (u) = 1 at 
u = w*-'^^ ~ (3.1,4.5) G Strip'-''"^^ this also clearly emerges from Figure [22] (in red). 

This is not sufficient to guarantee compatibility between the Ihs and rhs of ( 18. 7p . 
in fact also the part related to on the rhs must tend to one. This implies that 

F^^_-^{u) must have a critical point at u = ^ (3.1,4.5) G Strip which presence 
is clearly observable in Figure [20] (in red). 

Finally, the reader can see that the function F^{u) plotted in Figure [22] vanishes 
at w = u^~^^ ~ (2.8,4.7) G Strip (in blue), this is related to a critical point u^^^ ~ 
(2.8,6.7) G Strip "^^^^ (in red) observable in F^{u) plotted in Figure [211 Notice that the 
picture emerging from this quick inspection matches perfectly the complexes of critical 
points defined in Proposition 14.11 The interest for the classification of these complexes 
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or strings is not purely academic but it is actually an important preliminary step for 
the classification of the excited states of the theory. 





2 2.5 3 3.5 4 4.5 5 

Re(u) 



Figure 20. Details of F^.Ju): L = 4, g = 0.5, fiy = 0. 
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2 2.5 3 3.5 4 4.5 5 

Re(u) 



Figure 22. Details of F^{u): L = A, g = 0.5, iiy = 0. 



§ 9. Conclusions 

The correspondence between strings and gauge theories provides exactly soluble 
models with remarkable properties. These important systems can be studied using the 
powerful tools developed over the years in the integrability context for the study of 
models mainly relevant to condensed matter physics in low dimensions. 
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A main objective of this article is to describe preliminary numerical results con- 
cerning a new variant of the TBA equations for the AdS5/CFT4 correspondence in the 
attempt to reveal the analytic properties of their solutions. 

Studying the evolution of the ground state energy at different values of the scale 
and the chemical potential, we have checked that -as predicted by supersymmetry- the 
equations correctly lead to a vacuum state with zero energy as fXy — j- in at arbitrary 
scale L. 

It was then possible to explore the functions Fa^"^ in the complex plane of the 
rapidity, allowing us to study the zeroes and the singularities of the functions Ya{u) 
and 1 -|- Ya{u). We have seen that the basic Y-system relates points between adjacent 
strips of the complex plane and verified the presence of square root branch cuts exactly 
at the points predicted by [H] and summarized in Table [TJ These results confirm that 
the analytic properties of the AdSs / CFT4 thermodynamic Bethe Ansatz equations are 
very different from those of the relativistic integrable quantum field theories. The Y 
functions live on a complicated and -up to now- only superficially explored Riemann 
surface with an infinite number of square root branch points. 

Finally our set of programs should be easily adaptable to other non-linear integral 
equations, as for example those for the AdS4/CFT3 correspondence, linking the type IIA 
string theory on AdS4 x CP'^ to the Chern-Simons model with supersymmetry Af = 6 
in three dimensions [531 EH ES] . 
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§ A. The S-matrix elements 



TBA equations (OEH). 



(A.l) 



Sy,Q{u, Z) = SQ^y{z,u) 



( 



x{z - ^Q) - x{u) \ xjz + ^Q) 
x{z + ^Q) - x{u) J xiz-^^QY 
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(A.2) 



■ xiz-^Q)-x{u + ^^M) \ /xiz+^ 
x{z + ^Q)-x{u + j;M)j [x{z-j^Q) 

M-l 



X 



x{z-^^Q)-x{u-^M) \ t^^ z-u-^^{Q-M + 2j) 
x{z + ^^Q)-x{u-^M)) yz-u+^^{Q-M + 2j) 



(A.3) Sm{u) 



g 

9 



(A.4) 

Skm{u) = JJ^ JJ^ S'((i^+2-2fc)-(M-20)('") 



K M 



k=l 1=1 

'u - - M\\ (u-^{K + M)\ """(-f^^)-! I u - ^{\K - M\ + 2k)^ ^ 



n 



u+^K-M\ \u+^{K + M) J-l \uj^L{\K-M\+2k) 

\ 9 / \ 9 / k, — 1 \ 9 

The elements Sq/q are: 

(A.5) S^'q{z^ u) = {Sq>q{z - w))"^(Eq/q(z, 

where Sq'q is the improved dressing factor for the mirror bound states 

„ , . FT IT ( ~ ^(^+t(Q'+2-2fc))x(^+^(Q-2;)) ^ 
(A.6) ^Q'q{z, ^) = 11 11 . 1 (^Q'Qiz. u), 

k=ll=l \ x{z+^{Q'-2k))x{u+j{Q+2-2l)) J 

with ctq'q evaluated in the mirror kinematics. A precise analytic expression for the 
mirror improved dressing factor has been given in [43], and a more compact integral 
representation in |37]. The equivalence between these two results was proved in [H] . 
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